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ABSTRACT 
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computing  the  mean  power  pattern.  For  both  the  MM  and  PO 
based  simulations  a  plate  in  the  xy  plane  is  represented 
by  a  collection  of  small  triangles .  The  roughness  of  the 
surface  is  obtained  by  randomly  displacing  the  z 
coordinates  of  the  triangle  nodes.  The  main  focus  of  the 
study  was  to  look  at  the  effects  rough  surfaces  have  on 
the  off  principal  plane  regions. 
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I.  INTRODUCTION 


Since  the  invention  of  radar  in  World  War  II  many 
applications  have  been  found  for  its  use  in  both  the  military 
and  civilian  sectors.  Military  radar  applications  include  the 
detection  and  tracking  of  ships,  aircraft  and  ground  vehicles. 
Modern  radar  systems  make  it  virtually  impossible  for  ships 
and  aircraft  to  approach  targets  without  being  detected  at 
great  distances.  For  this  single  reason,  the  military  has  put 
great  emphasis  on  the  lowering  of  the  radar  cross  section 
(RCS)  of  military  vehicles.  The  word  stealth  has  been  used  by 
the  media  to  describe  specialized  military  vehicles 
specifically  designed  with  a  low  RCS.  In  principle,  a  stealth 
vehicle  can  approach  a  radar  without  being  detected,  and  thus 
has  an  increased  probability  of  completing  its  mission 
objectives  relative  to  a  conventional  vehicle. 

Several  factors  contribute  to  the  RCS  of  a  vehicle.  The 
shaping  of  its  surfaces  and  the  materials  used  in  its 
construction  have  the  most  significant  effect.  The  RCS  of  a 
rough  surface  can  be  significantly  different  from  that  of  a 
smooth  surface.  A  randomly  rough  surface  reduces  specular 
lobes  and  increases  diffuse  scattering.  Diffuse  scattering  is 
not  desirable  because  it  is  isotropic  and,  even  though  it  has 
a  very  low  level,  it  cannot  be  controlled  by  shaping.  Surface 
roughness  also  affects  the  behavior  of  surface  waves. 

The  sources  of  target  surface  imperfections  are 
unavoidable.  They  arise  from  tolerances  in  the  machining  and 
assembly  of  parts,  and  in  the  fabrication  of  materials 
themselves.  For  example,  the  wing  of  an  aircraft  typically 
consists  of  a  skin  material  laid  over  support  structures.  The 
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skin  will  have  some  deviation  from  the  theoretical  wing  shape. 
Furthermore,  if  the  skin  is  a  composite  material,  its 
dielectric  constant  and  loss  tangent  will  not  be  precisely  the 
same  throughout  the  material.  Thus  the  wing  is  a  rough, 
inhomogeneous,  anisotropic  scatterer. 

RCS  prediction  methods  can  be  divided  into  two 
categories:  1)  approximate  methods  such  as  physical  optics 
(PO)  and  2)  rigorous  methods  such  as  integral  equation 
techniques  and  the  method  of  moments  (MM) .  In  this  thesis 
both  PO  and  MM  are  used  to  calculate  the  RCS  of  plates  of 
various  sizes,  both  perfectly  flat  and  rough.  For  both  of 
these  methods  the  plate  is  represented  using  triangular 
patches.  When  MM  is  applied,  the  patches  are  the  subdomains 
used  in  the  current  expansion.  For  the  PO  solution,  closed- 
form  expressions  for  scattering  from  a  triangular  plate  are 
employed.  The  rough  surface  is  obtained  by  selecting  a  random 
number  at  each  node  of  the  plate.  The  random  deviations 
generated  in  this  manner  are  scaled  to  produce  a  random 
surface  of  a  specified  variance.  The  usefulness  of  PO  can  be 
accessed  by  a  comparison  of  the  average  RCS  level  with  that 
obtained  using  MM.  It  is  the  average  RCS  level  (or  diffuse 
scattering)  that  is  of  interest  in  specifying  tolerances  on 
the  manufacturing  and  asseinbly  processes. 

This  thesis  investigates  the  characteristics  of 
electromagnetic  scattering  from  surfaces  with  random  roughness 
and  its  impact  on  RCS.  The  first  objective  is  the  comparison 
of  the  approximate  and  rigorous  scattering  analysis  methods 
using  a  Monte  Carlo  simulation.  The  average  RCS  levels 
obtained  from  the  Monte  Carlo  simulations  are  then  compared  to 
a  closed-form  expression  for  the  mean  power  pattern  derived  by 
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Ruze  (Ruze,  1952) .  The  second  objective  is  to  investigate  the 
effects  of  surface  roughness  on  the  RCS  patterns;  in 
particular,  how  it  affects  the  specular,  traveling  wave,  and 
diffuse  RCS  components. 

Chapter  II  is  a  review  of  the  Rayleigh  Criterion  and  how 
rough  surfaces  affect  specular  and  diffuse  scattering.  Also 
discussed  is  a  method  for  the  calculation  of  the  mean  power 
pattern  using  a  statistical  approach. 

Chapter  III  reviews  PO  and  MM.  The  analysis  methods  are 
presented  and  the  required  equations  derived.  Computer  coding 
of  the  equations  in  MATLAB  is  also  discussed. 

Chapter  IV  is  an  evaluation  of  the  data  obtained  from  the 
computer  models,  primarily  highlighting  the  differences  in  RCS 
patterns  using  the  two  analysis  techniques.  A  summary 
consisting  of  the  mean  power  patten  difference  as  related  to 
the  roughness  of  the  surface  is  presented.  Chapter  V 
siammarizes  the  results  and  highlights  the  advantages  and 
disadvantages  of  each  approach  as  well  as  suggests  possible 
future  work. 


3 


n.  THEORETICAL  BACKGROUND 


In  this  chapter  the  concepts  and  equations  used  to 
describe  scattering  from  rough  surfaces  are  introduced. 

A.  RAYLEIGH  CRITERION 

One  of  the  first  researchers  to  study  the  effects  of 
scattering  from  rough  surfaces  was  Rayleigh  back  in  1877.  His 
approach  was  to  look  at  a  plane  wave  normally  incident  onto  a 
sinusoidal  surface.  [Ref.  1]  This  approach  led  to  the  basis 
for  determining  the  degree  of  roughness  of  a  surface,  now 
known  as  the  Rayleigh  criterion.  Figure  1  shows  a  rough 


Figure  1.  Rough  surface  diagram  used  to  determine  the  phase 
difference  between  two  parallel  scattered  rays.  After  Ref. 
[1]  . 
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surface  with  a  plane  wave  incident  from  angle  61  in  the  x-z 
plane.  The  phase  difference  (A4>)  between  the  two  parallel 
rays  scattered  into  the  x-z  plane  from  separate  points  can  be 
calculated  by 


A<l>=;t[(Zj-Z2)(cosej+cose2)+(x2-Xj)(sin0j-sine2)]  (2.1) 

where  the  heights  of  the  two  separate  scattering  points  are  Z/ 
and  Z;.  In  equation  2.1  the  variable  k  represents  the  modulus 
of  the  incident  and  scattered  wave  vector  (A=2nA)  .  For 
specular  scattering  6i=02  and  the  phase  difference  becomes 

A4>=2^(Zj-Zj)cos6j  .  (2.2) 

The  precise  way  in  which  the  two  scattered  rays  interact 
with  each  other  depends  on  the  value  of  Ad  relative  to  n.  If 
Ad  «  n,  then  the  two  scattered  waves  are  nearly  in  phase  and 
thus  interfere  constructively;  if  Ab  *  n,  then  the  two 
scattered  waves  interfere  destructively.  In  the  latter  case 
there  will  be  no  contribution  of  scattered  energy  in  the 
specular  direction.  [Ref.  1]  The  Rayleigh  criterion  is  a 
commonly  accepted  definition  of  a  rough  surface:  if  AcJ)  <  n 
then  the  surface  is  considered  smooth,  otherwise  it  is  rough. 
For  a  continuous  surface  with  many  scattering  points  the  root 
mean  square  (RMS)  deviation  5  is  used  in  place  of  ZyZ^.  Thus  the 
Rayleigh  criterion  can  be  expressed  as 


R  <ir/4 

a 


(2.3) 
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where  is  known  as  the  Rayleigh  parameter  and  is  given  by 


R  =kZcosQ.  . 

a  1 


The  Rayleigh  criterion  can  be  used  to  illustrate  the 
effects  of  a  rough  surface  on  an  incident  plane  wave.  If 
Zj=Z2,  then  the  phase  difference  becomes 


Figure!,  a)  Scattering  from  smooth  surface,  b)  scattering 
from  a  slightly  rough  surface,  and  c)  scattering  from  a 
rough  surface.  These  diagrams  show  how  the  degree  of 
surface  roughness  reduces  the  scattering  in  the 
specular  direction. 


A(()=A:(x2-Xj)(sm6j-sm62)  . 


(2.5) 


As  shown  in  Figure  2,  for  a  smooth  plate  of  very  large  extent, 
the  incident  wave  is  almost  completely  scattered  in  the 
specular  direction.  However,  if  Zj  is  not  equal  to  the 
surface  becomes  rough  and  less  energy  is  scattered  in  the 
specular  direction.  Equation  2.2  can  be  used  to  find  the 
phase  difference  for  the  given  surface  and  incident  angle. 
When  is  small  compared  to  n,  then  the  surface  can  be 
treated  as  though  it  were  smooth.  If  the  phase  difference  is 
no  longer  negligible,  a  significant  amount  of  the  incident 
energy  is  scattered  diffusely.  Finally,  if  the  phase 
variation  from  point-to-point  on  the  surface  becomes  large 
then  diffuse  scattering  dominates.  A  perfectly  diffuse 
surface  is  one  that  scatters  uniformly  in  all  directions  for 
every  incidence  angle. 

B.  STATISTICAL  APPROACH  TO  SOLID  SURFACE  SCATTERING 

Ruze  investigated  the  effects  of  random  surface  errors  on 
the  gain  of  reflector  antennas  using  a  statistical  approach. 
[Ref.  2]  It  is  possible  to  derive  a  closed-form  expression 
for  the  mean  power  pattern  from  the  radiation  integrals  for 
electromagnetic  scattering.  To  derive  the  expression,  some 
simplifying  assumptions  and  approximations  are  made.  One  is 
that  the  surface  bumps  can  be  represented  locally  by  a 
Gaussian  shape  as  shown  in  Figure  3.  The  width  of  the 
Gaussion  shape  is  characterized  by  a  correlation  interval,  c. 
It  is  the  average  distance  at  which  the  deviation  becomes 
uncorrelated.  An  easy  way  to  explain  large  correlation 
interval  is  to  imagine  a  slowly  varying  deviation  as  an 
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observer  travels  along  the  surface.  On  the  other  hand,  a 
small  correlation  interval  has  a  rapidly  varying  surface 
profile.  Figure  3  gives  examples  of  large  and  small 
correlation  intervals. 


S^<  variance 


Large  Correlation  Interval 

note  gradual  change  along  x  axis 


IkJS _ . 

r  V 

—  Gaussic 

1  V 

in  dent 

Small  Correlation  Interval 
note  rapid  changes  along  x  axis 

Figures.  Rough  surfaces  with  a)  large  correlation  interval 
and  b)  small  correlation  interval. 


A  closed-form  expression  for  the  average  RCS  can  be 
derived  from  the  mean  power  pattern 


<  p  >=<  E^-E*  > 


(2.6) 


where  <•>  denotes  expectation.  The  complete  derivation  for 
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the  mean  power  pattern  can  be  found  in  Jenn  [Ref.  3]  and 
therefore  only  the  highlights  are  presented  here.  Along  with 
the  gaussian  approximation  just  described,  a  second  assumption 
in  the  derivation  of  the  closed  form  equation  is  that,  as  long 
as  the  roughness  on  the  surface  is  sufficiently  small,  the  z 
component  of  the  surface  current  for  a  plate  lying  in  the  x-y 
plane  can  be  neglected.  Only  the  x  and  y  components  contribute 
to  the  scattered  field.  For  an  incident  wave  polarized 
transverse  magnetic  to  the  r-axis  (TM^) 


(2.7) 

- 

*  O 

where  g=-k»r',  k=k(itu+S>v-^tw) ,  and  r'=Stx'  andpy' .  The  quantities  u, 
V,  and  w  are  the  direction  cosines  of  the  observation  point 
and  are  defined  as:  u  =  sinGsincJ),  v  =  sin6sin(}),  w  =  cos6,  and 
(0,4))  gives  the  direction  to  the  observation  point.  rio 
represents  the  impedance  of  free  space. 

The  scattered  field  for  a  unit  incident  field  (£oe=l)  is 

E (2.8) 

where  represents  the  distance  from  the  local  surface  origin 
to  the  observation  point.  The  surface  devation  as  a  function 
of  position  is  5(r').  Now  using  equation  2.8  in  2.6 
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<  p  >=<  E  ‘E  >= 

S  S 


^COS0 


\  2 

ffj >ds 'ds.  (2.9 


The  error-free  power  pattern  (i.e.  pattern  of  a  prefectly  flat 
plate)  is 


icosd 


(2.10) 


Taking  the  expected  value  of  the  exponential  in  equation  2.9, 
expanding  the  integrand  in  a  Taylor  Series,  then  breaking  out 
the  first  term  gives 


<P>= 


^cosel ^  .4tS5sf  p  .  rr  Y' 


,  (2.11) 


where  and  6^  is  the  variance  of  the  deviations.  Switching 
to  polar  coordinates  and  extending  the  limits  of  the 
integration  to  infinity  for  reasons  stated  by  Jenn  [Ref.  3], 
equation  2.11  can  be  further  reduced  to 


<  p  >=  p  f  _ 


c^jt^sin^0 


(2.12) 


where  A  is  the  area  of  the  plate. 

The  RCS,  a,  can  now  be  found  in  terms  of  the  mean  power 
pattern  using  the  fact  that  the  incident  plane  wave  has  unit 
amplitude.  For  large 


o  =4jtrp <  P  > 


(2.13) 
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Equation  2.10  can  be  normalized  to  the  area  squared 


and  therefore  equation  2.13  becomes 


(2.14) 


4ltA  cos  6  n 

O  = - e  S  -r. 


+ - exp 

‘  A 


c^n^sin^d 


(2.15) 


The  second  term  in  equation  2.15  represents  the  average  or 
diffuse  RCS  level.  Because  the  variance  is  present  in  the 
second  term,  it  shows  the  importance  of  maintaining  tight 
surface  tolerances. 

A  similar  calculation  can  be  performed  for  a  two 
dimensional  array  of  discrete  isotropic  scatterers  that  are 
randomly  dislocated  from  the  z=0  plane.  The  mean  power 
pattern  consists  of  two  terms  similar  to  equation  2.15.  The 
first  is  the  error-free  pattern,  which  depends  on  the  array 
factor  of  the  scatterers.  The  second  term  is 


4iiA^  cos^0(4fcos?05^) 
MN 


where  MN  is  the  total  number 
array.  Equation  2.16  represents 


of  scatterers  making  up  the 
the  average  RCS  level.  It  is 
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interesting  to  note  that  if  the  discrete  array  has  a 
rectangular  grid  spaced  dx  by  dy,  its  area  is  MNdxxfy.  Equating 
the  second  term  of  equation  2.15  to  equation  2.16  gives 


dxdy=%c^  . 


(2.17) 


Thus  an  equivalent  correlation  interval  for  the  discrete  array 
can  be  computed. 
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m.  METHOD  AND  FORMULATIONS 


The  equations  required  for  the  PO  and  MM  Monte  Carlo 
simulations  are  developed  in  this  chapter.  Also  introduced  is 
a  brief  description  of  the  computer  code  used  in  the 
simulation. 

A.  PHYSICAL  OPTICS  (PO) 

Probably  the  most  commonly  used  prediction  method  for  the 
study  of  wave  scattering  from  rough  surfaces  is  the  physical 
optics  approximation.  [Ref.  1]  Physical  optics  provides  an 
estimate  of  the  surface  current  on  the  illuminated  part  of  the 
target.  The  shadowed  area  is  assumed  to  have  a  surface 
current  of  zero.  One  limiting  factor  that  should  be  mentioned 
is  that  the  use  of  physical  optics  will  not  show  the  effects 
of  traveling  waves  that  normally  appear  along  the  edges  of 
surfaces.  Surfaces  that  are  terminated  with  edges,  as  in  the 
case  of  a  rectangular  plate,  have  traveling  waves  reflected 
back  towards  the  source.  Thus  for  monostatic  RCS 
calculations,  PO  can  give  erroneous  results  at  angles  near 
surface  grazing. 

The  derivation  of  the  scattered  fields  for  physical 
optics  begins  with  the  geometrical  optics  approximation  of  the 
surface  current  on  the  illuminated  portion  of  the  body.  The 
surface  current  is 


j  ]2fi  X  ,  for  illuminated  part  of  the  body,  (3.1) 

'  [  0,  for  shadowed  part  of  the  body. 

In  equation  3.1,  H,-  is  the  incident  magnetic  field  intensity 
at  the  surface  and  ti  is  the  surface  normal  vector.  Now 
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consider  a  flat  plate  of  dimensions  a  and  b  as  shown  in  Figure 
4.  An  incident  wave  has  arbitrary  polarization  defined  by  the 
constants  and  Eo^,  and  is  incident  from  a  direction  given 
by  the  angles  (0,4)): 

,3.2) 


The  magnetic  field  intensity  can  be  written  as 


Figure  4.  Plane  wave  scattering  by  a  flat  plate. 
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provided  that  the  wave  is  propagating  towards  the  origin.  Now 
substituting  the  magnetic  field  intensity  into  equation  3.1  it 
is  possible  to  approximate  the  surface  current  of  the 
rectangular  plate 

X  .  (3.4) 


For  a  flat  plate  in  the  x  y  plane,  equation  3.4  becomes 


_  A  A  gA? 

J«2fx(0£^^-4,£^^) -  (3.5) 

L  . 

where  g=x'u+y'v+z'w. 

Finally,  the  monostatic  scattered  field  can  be  determined 
from  the  surface  current 


E^J^ff  j 
*  jJa  ' 


Anr^'' 


(3.6! 


Combining  equations  3.5  and  3.6  gives 


E=JL.e  -$£  ]) f  f  e^*»yds '  (3.7) 

where  A  is  the  total  area  of  the  plate  and  g=h  for  the 
monostatic  case. 

A  rough  plate  can  be  represented  by  a  collection  of 
triangular  patches.  The  edges  of  adjacent  triangles  touch,  but 


17 


the  normals  point  in  different  directions  as  evident  in  Figure 
5.  Thus  the  total  integral  in  3.7  is  a  sum  of  integrals,  each 
over  a  single  triangle 


(3.8) 


where  W,  represents  the  total  number  of  triangles  on  the  plate 
and  Ij  is  the  ith  integration  result.  The  closed-form 
expressions  for  /,■  are  derived  by  Moreira  and  Prata.  [Ref.  4] 
In  general. 


[Dp-D;i  ‘ 


D  D 

P  9 


(3.9) 


where  is  the  area  of  the  triangular  cell.  Figure  6  shows 
the  integration  cell  geometry  used  in  equation  3.9.  For 
plane  wave  incidence 

C=Aij-Ai3^0,  C^=A,.rA,-j=0.  Co=A^j,  F,=0, 

D=k{B,rB,s),  Do^kB,,,  F,=0. 

A^„  and  Bj„  represent  the  amplitude  and  phase  of  the  integrand 
at  node  m  of  triangle  /,  respectively,  as  shown  in  Figure  6. 
Equation  3.9  can  be  reduced  to 


D  (P  -D) 

P  9  P  ] 


D  (P  -D)\  D  D 

q  \  n  t 


(3.10) 


Equation  3.10  has  some  removable  singularities  that  cause 
problems  numerically  when  |T>p|,  ID^l  and  {D^-Dpl  are  small. 
These  singularities  are  handled  as  special  cases  and 
expressions  are  derived  by  expanding  the  integrand  in  a  Taylor 
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Figures.  Triangular  patch  model  of  a  rough  scattering 
surface . 
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series  and  analytically  integrating  the  resultant  series  term 
by  term.  The  following  special  cases  result  [Ref.  4]: 

1.  If  \Dp\<Lt  and  \D^\>L^  (no  restrictions  on  ID^-Dpl) 


Figure  6.  Integration  cell  geometry  in  local  coordinates 


2.  If  \Dp\<Lt.  and  |f>^|<Xt  (no  restrictions  on  \Dg-Dp\) 


B-O  m-O  («+w+2)! 


(3.12: 


3.  If  IDpULt  and  |D^|<Xt  (no  restrictions  on  \D^-Dp\) 


'qf  '"t 


n-O  «!  \  «  +  U 


G(«+1,-Z)p 


(3.13) 


4.  If  and  \D^UL,  and  |D^-DJ<i 


-Cfiitt,  D  )tC,- 


jD^n-0  ft 


(3.14) 


The  function  G  is  defined  as 


X 

G(n,a)=fs”e^"ds 


(3.15) 


and  can  be  evaluated  recursively  using  the  relationship 


G(n,a)= - i ,  fortikl  , 


(3.16) 


G(0,a)=i^^  . 

/« 


(3.17) 


As  noted  by  Moreira,  equation  3.16  is  numerically  efficient 
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and/  accurate  results  are  obtained  in  most  cases  by  setting 
Lt=.05  and  using  n,m^  2.  However  for  the  rough  plate  it  was 
found  that  Lt=.001  was  required  to  avoid  spikes  in  the  RCS 
pattern.  Finally,  the  RCS  of  the  plate  can  be  calculated 
from 


o  =4nrp  (3.18) 

and  by  using  equation  3.7  and  3.8  to  determine  the  scattered 
field. 

B.  COMPUTER  SIMULATION  USING  PHYSICAL  OPTICS 

The  rough  surfaces  used  for  the  computer  calculations  in 
this  thesis  were  created  using  MBUILD.  It  is  written  in 
FORTRAN,  and  can  generate  basic  surfaces  and  then  approximate 
them  with  triangular  facets.  The  output  is  a  data  file 
consisting  of  the  x,  y,  and  z  node  locations  of  each  facet. 
The  z  node  is  then  randomly  displaced  to  generate  a  rough 
surface  of  specified  standard  deviation.  For  the  physical 
optics  approximation,  this  data  file  is  then  broken  up  into 
various  MATLAB  files  consisting  of  individual  x-y-z  points, 
edges,  vertices,  and  faces. 

The  physical  optics  approximation  calculation  was  done  in 
MATLAB.  The  code  is  designed  to  present  the  user  a  three 
dimensional  plot  of  the  RCS  pattern  in  direction  cosine  space. 
Equations  3.10  through  3.14  are  functions  of  the  local 
coordinate  angles  and,  therefore,  to  use  these  equations,  a 
transformation  must  be  performed  for  each  triangle  as  shown  in 
the  Figure  7.  The  transformation  involves  rotating  about  the 
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z  axis  to  the  plane  defined  by  f  and  ft,  and  then  rotating 
through  the  angle  between  f  and  ti  in  that  plane. 


Figure?.  Translation  of  coordinates  for  a  plate  triangle. 


C.  METHOD  OF  MOMENTS  (MM) 

MM  is  a  numerical  technique  that  is  very  powerful  because 
it  can  be  applied  to  arbitrary  bodies.  Basically  MM  reduces 
an  integral  equation  to  a  matrix  problem.  The  matrix 
dimensions  are  dependent  on  the  size  of  the  scattering  body  in 
terms  of  the  incident  wavelength.  For  this  reason  the 
application  of  MM  is  limited  to  small  and  moderate  target 
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sizes . 


The  integral  equation  to  be  solved  is  the  E-field 
integral  equation  (EFIE) 

I f  F')-— [V'1,(F')]V'G(F,  .  (3.19) 

^  j  j  s  •  Qe 

In  this  chapter,  S  denotes  the  surface  of  integration. 

The  MM  solution  begins  by  representing  the  surface  current 
by  a  series  with  unknown  expansion  coefficients,  I„  : 

<3-20) 

#1*1 

J  are  the  basis  function,  and  both  I„  and  J„  can  be  complex. 
An  important  consideration  in  the  MM  solution  is  the  selection 
of  the  basis  functions.  One  should  chooses  a  basis  function 
that  has  the  behavior  expected  of  the  current  and  minimizes 
computational  effort.  These  functions  are  typically  pulses, 
triangles,  sinusoids  and  5  functions. 

The  next  step  is  to  solve  for  the  series  coefficients. 
Define  a  set  of  weighting  functions  W„,  /w=l,2,...N.  It  is 
convenient  to  choose  the  complex  conjugate  of  the  basis 
function  being  used;  this  choice  is  referred  to  as  Galerkin's 
method.  [Ref.  3]  The  EFIE  is  multiplied  by  each  weighting 
function  and  integrated  over  the  domain  resulting  in  N 
equations  of  form 
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/  is  a  vector  containing  the  unknown  expansion  coefficients. 
Z is  the  impedance  matrix  and  V  is  the  excitation  vector.  By 
inverting  the  matrix  Z,  /  can  be  solved  for: 

/=Z‘^F.  (3.26) 


Having  evaluated  the  equation  3.20  is  now  defined  and  can 
be  used  in  the  radiation  integral  (equation  3.6) . 


A  measurement  matrix,  R  can  be  defined  with  elements 


(3.27) 


where  E,  is  a  unit  radiated  plane  wave.  [Ref.  3] 
polarized  wave 


r 


Using  a  0 


[3. 28) 


Combining  equations  3.27  and  3.28 


rK 1 1  [e-J 


(3.29) 


where  the  superscript  0  is  added  to  denote  the  polarization  of 
E,.  Note  that  u,  v,  and  w  are  the  direction  cosines  for  the 
observation  direction.  The  same  can  be  done  for  a  (p  polarized 
plane  wave.  The  electric  field  components  are 


-Jkr\  ^  P  a 


(3.30) 


and  for  4) 


Er-^e*‘ZRX 


(3.31) 


Referring  to  Figure  5,  the  rectangular  plate  is  broken  up 
into  a  collection  of  triangular  faces-  The  surface  current  is 
represented  by  equation  3.20  where  r  is  a  position  vector  to 
a  source  point.  Figure  8  shows  the  patch  geometry  for  two 
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triangles  with  one  common  edge.  The  niomber  of  basis  functions 
used  is  dependent  on  the  number  of  edges  shared  with  other 
triangles.  If  a  patch  has  two  shared  edges,  then  it  is 
spanned  by  two  basis  functions;  a  patch  with  three  shared 
edges  would  be  spanned  by  three  basis  functions.  In  Figure  8, 


Figures.  Triangular  patch  geometry. 


edge  n  is  shared  by  two  triangular  patches,  In'"  and  T„’  with 
areas  and  A„'.  The  shared  edge  has  a  length  of  L„.  The 
vector  basis  function  is  written  as 
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F  in  T 


(3.32) 


r  in  T' 


Some  characteristics  of  these  functions  are:  1)  for  triangular 
edges  that  represent  body  edges,  there  is  no  normal  component 
of  the  surface  current,  2)  the  current  normal  to  an  internal 
edge  n  is  constant  and  continuous  across  the  edge,  and  3)  the 
divergence  of  J„  exists  and  has  the  expression 


F  in  T: 


(3.33) 


r  in  T, 


The  MM  impedance  elements  can  now  be  written  as 


(3.34) 


where 


(3.35) 


=_ll 

mn  . 


4llJ(x>e 


rr  - 

//vV.(r-')i— 


(3.36) 
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where  R*'=\r‘^-r'\.  The  excitation  vector  elements  are 

m  I  m  • 


V  =z  (£*-pr+£^'*pr) 

IN  in'*  m  m  r  m  / 


(3.37) 


where  the  incident  field  is  evaluated  at  the  centroid 


(3.38) 

Using  the  expression  for  the  impedance  elements  in  3.34, 
the  integrals  must  be  evaluated  numerically.  This  should  be 
accomplished  by  transforming  to  a  set  of  localized  coordinates 
on  each  patch  and  defining  a  set  of  normalized  coordinates  as 
in  the  case  of  PO.  Singularity  correction  of  the  self  term 
(m-n) of  the  impedance  matrix  is  generally  required. 

D.  COMPUTER  SIMULATION  USING  MM 

The  MM  code  used  to  simulate  the  varying  rough  surface  was 
a  code  written  by  Sandia  National  Laboratories  called  PATCH. 
[Ref.  5]  The  output  from  the  code  provides  the  user  with  the 
currents  coefficients  and  radar  cross  sections.  PATCH  uses 
the  geometry  file  generated  by  MBUILD,  which  was  previously 
discussed.  Minor  modifications  were  made  to  PATCH  to 
calculate  the  RCS  pattern  of  the  rough  surfaces  in  direction 
cosine  space  for  comparison  with  the  PO  results. 
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IV.  SUMMARY  OF  RESULTS 


A.  ROUGH  SURFACE  PLATE  MODELS 

The  surfaces  considered  in  all  the  simulations  are  3X  by 
3X  square  plates.  This  was  the  largest  plate  possible  given 
computer  limitations.  A  large  plate  is  desirable  so  that  the 
principle  plane  sidelobes  fall  off  rapidly,  thus  allowing  the 
diffuse  RCS  of  the  off  principal  planes  to  dominate.  Another 
limitation  was  the  computer  run  time. 

The  FORTRAN  code  MBUILD  was  used  to  construct  a  plate 
comprised  of  1800  triangles  as  shown  in  Figure  9.  The  plate 
node  coordinates  were  then  altered  to  produce  the  rough 
surfaces  shown  in  Figure  10  and  Figure  11.  Figure  10  has  a 
surface  deviation  uniformly  distributed  from  0  to  0 .  lA  with  a 
mean  of  0.05X  and.  Figure  11  has  deviation  varying  from  0  to 
0.2A  with  a  mean  of  O.lX.  The  reason  for  the  O.IX  and  0.2X 
deviations  was  to  create  a  significant  effect  on  the  RCS 
pattern  in  the  off  principal  plane  regions.  The  manipulation 
of  the  z  coordinates  was  done  using  a  FORTRAN  program  ZCHN 
which  generates  a  uniform  random  number  from  0  to  1  and  then 
is  scaled  based  on  the  desired  maximum  amplitude  deviation. 
Thus  three  plates  are  used  in  the  simulation  models  in  order 
to  study  their  effects  on  RCS  patterns. 
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Figure  9.  Flat  plate  3.  OX  by  3.  OX. 
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Figure  10.  Rough  plate  with  O.lX  maximum  surface  deviation. 


Figure  11.  Rough  plate  with  0.2X  maximum  surface  deviation. 
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B.  METHOD  OF  MOMENTS  RESULTS 

Method  of  moments  is  a  rigorous  solution  of  the  scattering 
equations  and  therefore  includes  the  effects  of  traveling 
waves.  The  software  used  to  calculate  the  RCS  patterns  was 
PATCH.  Figure  12  shows  the  RCS  pattern  for  a  flat  3A  by  3A 
plate.  The  traveling  waves  cause  the  high  RCS  levels  present 
around  the  edges  of  the  pattern.  The  main  lobe  and  principal 
plane  side  lobes  are  quite  noticeable,  but  the  off  principal 
plane  lobes  are  not  as  well  defined  because  of  the  traveling 
wave.  When  the  surface  roughness  is  increased  to  O.lA,  as  in 
Figure  13,  the  principal  plane  lobes  become  less 
distinguishable  because  diffuse  scattering  raises  the  average 
level.  Increasing  the  roughness  to  0.2X  ,  as  shown  in  Figure 
14,  results  in  a  relatively  uniform  RCS  level  everywhere 
except  the  main  lobe.  Even  the  traveling  wave  effects  have 
been  reduced  by  the  rough  surface.  In  Figure  15,  a  cut  is 
shown  for  u=-0.27.  It  illustrates  the  effect  of  the  surface 
roughing:  the  principal  plane  lobes  are  decreased  and  the  off 
principal  plane  lobes  are  increased. 
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Figure  12.  RCS  of  a  flat  plate  computed  using  method  of 
moments  (3A  by  3A) . 


rough  plate  co: 


C.  PHYSICAL  OPTICS  RESULTS 

Using  the  MATLAB  code  POCOS,  the  RCS  of  3X  by  3A  plates 
where  plotted  in  direction  cosine  space  for  comparison  with 
the  MM  data.  Figure  16  shows  a  RCS  pattern  for  a  flat  plate 
using  the  PO  approximation.  Because  PO  does  not  include 
traveling  waves,  the  RCS  levels  around  the  edge  of  the 
pattern  are  noticeably  lower  than  those  of  MM.  The 
periodicity  of  the  lobes  is  apparent  when  the  RCS  is  plotted 
in  direction  cosine  space.  Figures  17  and  18  show  the  RCS 
patterns  for  the  rough  plates  =  O.IA  and  0.2A)  .  The 
effects  of  random  errors  are  more  noticeable  in  this  case  than 
they  are  for  MM  because  of  the  absence  of  traveling  waves.  As 
the  surface  roughness  increases,  the  off  principal  plane  lobes 
increase,  and,  if  the  roughness  increases  without  bound,  the 
scattering  from  the  plate  would  become  completely  diffuse.  A 
comparison  of  the  rough  surface  patterns  are  shown  in  Figure 
19.  A  cut  of  the  RCS  pattern  has  been  taken  at  u=-0.27. 
Figure  19  shows  the  effect  of  the  surface  roughing.  The  off 
principal  plane  levels  in  some  areas  increase  by  as  much  as 
35dB,  while  principal  plane  levels  decrease  only  slightly. 
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Figure  16.  RCS  of  a  flat  plate  computed  using  physical  optics 
(3X  by  3X)  . 
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D.  COMPARISON  OF  THE  METHODS 

The  closed-form  approximation,  equation  2.15,  can  be  used 
to  calculate  the  average  RCS  level  based  on  the  area, 
variance,  and  correlation  interval  of  the  plate  surface.  The 
second  term  in  equation  2.15  is  plotted  in  Figures  20  and  21. 
The  curves  for  several  correlation  intervals  are  shown.  The 
MM  and  PO  RCS  curves  are  averaged  over  the  four  quadrants. 
Furthermore  values  of  |u|<0.2  and  |v|<0.2  where  discarded  when 
averaging,  because  the  RCS  in  these  regions  is  dominated  by 
the  specular  RCS  (ie.,  first  term  in  equation  2.15)  .  Surface 
roughness  is  not  as  noticeable  in  the  MM  plot  because  the 
surface  wave  scattering  dominates  at  wide  angles,  even  off  the 
principal  planes.  The  PO  patterns  illustrate  that  the  diffuse 
scattering  does  dominate  over  the  specular  scattering  at  wide 
angles  (0^40*)  . 

The  plate  nodes  are  located  every  O.IA  in  the  x  and  y 
directions,  and  the  roughness  is  achieved  by  randomly  choosing 
the  z  coordinate  at  each  node.  Therefore,  a  first  order 
estimate  of  the  correlation  interval  might  be  c=0.lX.  This 
value  gives  curves  too  high  relative  to  PO,  which  is  the  best 
indicator  of  the  diffuse  level  (because  it  has  no  traveling 
waves) .  From  the  curves  shown  in  the  figures,  it  appears  that 
the  correlation  interval  is  actually  less,  perhaps  0.07X. 
This  is  close  to  the  result  obtained  equating  the  second  terms 
of  equations  2.15  and  2.17.  From  equation  2.17 

_2_(o.ony 

or  c~0.06X  . 
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Note  that  the  error  terms  represent  the  average  levels  due  to 
surface  roughness.  The  RCS  has  a  Rayleigh  probability  density 
function  and  it  is  expected  that  ninety-eight  percent  of  the 
sidelobes  will  be  less  than  this  average  value  plus  6dB. 
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Figure  20.  Comparison  of  the  three  prediction  methods  (MM,  PO, 
and  closed-form)  for  a  3X  by  3K  plate  and  O.lA. 
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Figure  21.  Comparison  of  the  three  prediction  methods  (MM,  PO, 
and  closed-form)  for  a  3X  by  3X  plate  and  5max=  0.2X. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


For  both  MM  and  PO,  as  the  surface  roughness  increased, 
so  did  the  RCS  levels  in  the  off  principal  planes.  Because 
MM  is  rigorous  and  takes  into  account  traveling  waves,  it 
represents  the  most  accurate  modeling  method.  It  yields 
higher  RCS  levels  in  the  off  principal  plane  regions  than  PO. 
However  MM  has  several  limitations  that  should  be  considered 
before  using  it.  Because  MM  is  so  computationally  intensive, 
very  fast  computers  with  large  amounts  of  memory  are  required, 
thus  limiting  the  size  of  object  that  can  be  modeled.  Another 
factor  is  the  amount  of  time  one  is  willing  to  wait  for 
results.  The  plates  modeled  in  this  thesis  took  approximately 
eight  days  of  computer  run  time  per  plate  on  a  SparclO  with  64 
Megabytes  of  memory.  MM  does  provide  the  most  accurate 
results  but,  on  the  other  hand,  the  PO  calculations  take 
approximately  one  day.  If  an  approximate  closed-form 
expression  could  be  found  for  the  traveling  wave  that  could  be 
added  to  the  PO  RCS  patterns,  the  PO  approximation  would  be 
the  best  choice,  yielding  results  similar  to  those  of  MM. 

Ruze's  closed-form  expression  can  be  used  to  predict  the 
RCS  in  regions  where  diffuse  scattering  dominates.  For  the  3X 
square  plate  there  is  really  no  angular  sector  where  diffuse 
scattering  dominates,  and  therefore  the  agreement  with  MM  is 
not  good.  However,  PO  does  not  have  a  traveling  wave  lobe  and 
therefore  diffuse  scattering  does  dominate  at  wide  angles  off 
the  principal  planes.  In  this  case  the  closed-form  results 
agree  with  PO  if  the  correct  correlation  interval  is  chosen. 

Future  work  should  focus  on  finding  an  accurate  closed- 
form  expression  for  the  traveling  wave  component  of  RCS  that 
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can  be  added  directly  to  the  PO  RCS.  Thus  the  total  RCS  can 
be  expressed  approximately  by 


That  is,  the  total  RCS  is  a  simple  summation  of  specular  (PO) , 
diffuse  (DIF),  and  traveling  wave  (TW)  terms.  Equation  5.1 
assiimes  that  the  sum  of  the  fields  squared  is  approximately 
the  sum  of  the  powers  and  this  is  true  if  only  one  of  these 
terms  dominates  in  any  particular  direction. 


APPENDIX.  COMPUTER  CODE 


This  appendix  contains  the  code  developed  to  calculate  RCS 
using  the  physical  optics  approximation.  Code  JPOCOS 
calculates  the  RCS  pattern  of  a  surface  of  random  roughness. 
In  present  form,  the  code  calculates  the  scattered  RCS  levels 
as  a  function  of  u  and  v  in  increments  of  .022  or  2/90.  This 
creates  a  data  string  of  8281  points  which  can  be  manipulated 
to  create  a  matrix  of  91  by  91  data  points.  This  data  is  then 
plotted  versus  u  and  v  giving  the  three  dimensional  patterns 
presented  in  Chapter  V. 
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%  program  POCOS.m 

%  latest  version  7/3 1/95  »>  Uses  geometry  data  from  build  «< 

%  PO  scattering  from  triangles 

%  ***  has  global  to  local  coordinate  transformations  *** 

%****  vertices  must  be  entered  in  clockwise  order  **** 
clear  all 

%  Dlumination  flag:  =0  external  face  only 
iflag=l; 

%  Geometry  flag  for  special  test  cases 
Lt=.0001;  %  Taylor  series  region 
Nt=2;  %  Number  of  terms  in  Taylor  series 
Rs=.0;  %  Surface  resistivity  (normalized) 
wave=l; 
blc=2*pi/wave; 
rad=pi/180; 

%  Incident  wave  polarization 
Et=l+j*0;  %TM-z 
Ep=0+j*0;  %TE-z 

%  Wave  amplitude  at  all  vertices 
Co=l.; 

%  Load  data  arrays  generated  by  build 
load  xpts.m 
load  ypts.m 
load  zpts.m 
nverts=length(xpts); 
load  endl.m 
load  end2.m 
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nedges=length(end  1 ); 
load  nodel.m 
load  node2.m 
load  nodes  .m 
nfaces=length(node3); 

%  Load  vind  array  which  gives  the  three  nodes  for  each  triangle 
for  i=l::nfaces 

pts=[nodel(i)  node2(i)  node3(i)]; 
vind(i,:)=pts; 
end 

ntria=nfaces; 

nvert=nverts; 

x=xpts; 

y=ypts; 

z=zpts; 

%  Define  position  vectors  to  vertices 
for  i=l  invert 
r(i,:)=[x(i)  y(i)  z(i)]; 
end 

%  Get  edge  vectors  and  normals  from  edge  cross  products 
for  i=l;ntria 

A=r(vind(i,2),  :)-r(vind(i,  1), :); 

B=r(vind(i,3),  :)-r(vind(i,  1 
C=r(vind(i,3),  :)-r(vind(i,2), :); 

N(i,  :)=cross(A,B); 

%  Edge  lengths  for  triangle  "i" 
d(i,  l)=norm(A); 
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d(i,2)=norm(B); 

d(i,3)=norm(C); 

ss=.5*sum(d(i,:)); 

Area(i)=sqrt(ss*(ss-d(i,l))*(ss-d(i,2))*(ss-d(i,3))); 
Nn=norm(N(i, ; )); 

N(i,:)=N(i,;)/Nn; 
t(i,:)=cross([0  0  l],N(i,:)); 
beta(i)=acos(N(i,3)); 
alpha(i)=atan2(t(i,2),t(i,  1)+1  e-6)-pi/2; 
end 

%  Pattern  loop 
it=0; 

foru=-l;2/90:l 
for  v=-l:2/90;l 
it=it+l; 

Sth(it)=-60;  Sph(it)=-60; 
if(u^2+v^2)<=l 

phi=atan2(v,u+ 1  e-8); 

cp=cos(phi);sp=sin(phi); 

theta=asin(sqrt(u^2+v'^2)); 

w=cos(theta); 

st=sin(theta);ct=cos(theta); 

D0=[u  V  w]; 

uu=ct*cp;  w=ct*sp;  ww=-st; 

%  Spherical  coordinate  system  radial  unit  vector 
R=[u  V  w]; 
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%  Incident  field  in  global  Cartesian  coordinates 
e0(  1  )=uu*Et-sp  *Ep; 
eO(2)=w*Et+cp*Ep; 
eO(3)=ww*Et; 

%  Begin  loop  over  triangles 
sumt=0; 
sunip=0; 
for  m=l:ntria 

%  Test  to  see  if  fi'ont  face  is  illuminated 
ndotk=N(m,  :)*R'; 

%  if  ndotk  >=0 

%  "local"  direction  cosines 

ca=cos(alpha(m));  sa=sin(alpha(m));  cb=cos(beta(m));  sb=sin(beta(m)); 
Tl=[ca  sa  0;  -sa  ca  0;  0  0  1];  T2=[cb  0  -sb;  0  1  0;  sb  0  cb]; 
D1=T1*D0'; 

D2=T2*D1; 

u2=D2(l);  v2=D2(2);  w2=D2(3); 

%  Find  spherical  angles  in  local  coordinates 

st2=sqrt(u2'^2+v2'^2)*sign(w2);  %  this  makes  a  difference 

ct2=sqrt(l-st2^2); 

phi2=atan2(v2,u2+ 1  e- 1 0); 

th2=acos(ct2); 

cp2=cos(phi2);  sp2=sin(phi2); 

%  Phase  at  the  three  vertices  of  triangle  m;  monostatic  RCS  needs  "2" 
Dp=2*bk*((x(vind(m,  l))-x(\dnd(m,3)))*u+. . . 

(y(vind(m,  l))-y(vind(m,3)))*v+... 

(z(vind(m,  l))-z(vind(m,3)))*  w); 
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Dq=2*bk*((x(vind(m,2))-x(vind(m,3)))*u+. . . 
(y(vind(m,2))-y(vind(m,3)))*v+. . . 
(z(vind(m,2))-z(vind(m,3)))*w); 

Do=2*bk*(x(vind(m,3))*u+y(vind(m,3))*v+z(vind(m,3))*w); 
%  Incident  field  in  local  Cartesian  coordinates  (stored  in  e2) 
el=Tl*eO'; 
e2=T2*el; 

%  Incident  field  in  local  spherical  coordinates 
Et2=e2(l)*ct2*cp2+e2(2)*ct2*sp2-e2(3)*st2; 
Ep2=-e2(l)*sp2+e2(2)*cp2; 

%  Reflection  coefficients  (Rs  is  normalized  to  etaO) 
perp=-l/(2*Rs*ct2+l);  %local  TE  polarization 
para=0;  %local  TM  polarization 

if  (2*Rs+ct2)~=0  para=-ct2/(2*Rs+ct2);  end 
%  Surface  current  components  in  local  Cartesian  coordinates 
Jx2=ct2  *  (-Et2  *  cp2  *  para+Ep2  *  sp2  *  perp); 

Jy2=ct2  *  (-Et2  *  sp2  *para-Ep2  *cp2  *perp); 

%  Area  integral  for  general  case 
DD=Dq-Dp; 
expDo=exp(j  *Do); 
expDp=exp(j  *Dp); 
expDq=exp(j  *Dq); 

%icase(i,m)=-l; 

%  Special  case  1 

if  abs(Dp)  <  Lt  &  abs(Dq)  >=  Lt 
%  icase(i,m)=l; 
sic=0.; 
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for  n=0:Nt 

sic=sic+(j  *Dp)^fact(n)*(-Co/(n+ 1  )+expDq*(Co*gf(n,-Dq))); 
end 

Ic=sic*2*Area(m)*expDo^/Dq; 

%  Special  case  2 

elseif  abs(Dp)  <  Lt  &  abs(Dq)  <  Lt 
%  icase(i,m)=2; 
sic=0.; 
for  n=0:Nt 
for  nn=0;Nt 

sic=sic+(j*Dp)^n*(j*Dq)'^nn/fact(nn+n+2)*Co; 

end 

end 

Ic=sic*2*Area(m)*expDo; 

%  Special  case  3 

elseif  abs(Dp)  >=  Lt  &  abs(Dq)  <  Lt 
%  icase(i,m)=3; 
sic=0.; 
for  n=0:Nt 

sic=sic+(j  *Dq)'^n/fact(n)*Co*gf(n+l  ,-Dp)/ (n+1 ); 
end 

Ic=sic*2*Area(m)*expDo*expDp; 

%  Special  case  4 

elseif  abs(Dp)  >=  Lt  &  abs(Dq)  >=  Lt  &  abs(DD)  <  Lt 
%  icase(i,m)=4; 
sic=0.; 
for  n=0:Nt 
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sic=sic+0*DD)^n/fact(n)*(-Co*gf(n,Dq)+expDq*Co/(n+l)); 

end 

Ic=sic*2*Area(m)*expDo/j/Dq; 

else 

%  icase(i,m)=0; 

Ic=2*Area(m)*expDo*(expDp*Co/Dp/DD-expDq*Co/Dq/DD-Co/Dp/Dq); 
end  %  end  of  special  cases  test 
%  Scattered  field  components  for  triangle  m  in  local  coordinates 
Es2(l)=Jx2*Ic;  Es2(2)=Jy2*Ic;  Es2(3)=0; 

%  Transform  back  to  global  coordinates,  then  sum  field 
Esl=T2’*Es2'; 

EsO=Tl'*Esl; 

Ets=uu*EsO(l)+w*EsO(2)+ww*EsO(3); 

Eps=-sp*EsO(l)+cp*EsO(2); 

%  Sum  over  all  triangles  to  get  the  total  field 
sumt=sumt+Ets; 
sump=sump+Eps; 

%  end  %  end  of  illumination  test 
end  %  end  of  triangle  loop 

Sth(it)=10*logl0(4*pi*abs(sumt)^2+l.e-10); 
Sph(it)=10*logl0(4*pi*abs(sump)^2+l.e-10); 

end  %  end  of  if  statement 
end  %  end  of  pattern  loop  V 

end  %  end  of  pattern  loop  U 

Sth(:)=max(Sth(:),-60); 
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Sph(;)=max(Sph(:),-60); 


Sth21=Sth'; 

Sph21=Sph'; 

save  Sthll.m  Sth21  -ascii 
save  Sphll.m  Sph21  -ascii 
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